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We revisit the sum-rule extraction of the decay constants of the D, D s , B, and B s mesons from the 
two-point correlator of heavy-light pseudoscalar currents. We use the operator product expansion of 
this correlator expressed in terms of the MS heavy-quark mass, for which the perturbative expansion 
exhibits a reasonable convergence. Our main emphasis is laid on the control over the uncertainties 
in the decay constants, related both to the input QCD parameters and to the limited accuracy 
of the method of sum-rules. This becomes possible due to the application of our procedure of 
extracting hadron observables that involves as novel feature dual thresholds depending on the Borel 
CNj , parameter. For charmed mesons, we find the decay constants fr> = (206. 2±7.3(ope) ±5.1( syst )) MeV 

and fo s — (245.3±15.7(ope) ±4.5( syst )) MeV. For beauty mesons, the decay constants turn out to be 
extremely sensitive to the precise value of nib (rn~b)- By requiring our sum-rule estimate to match the 
average of the lattice results for fs, a very accurate value = 4.245±0.025 GeV is extracted, 

leading to f B = (193.4 ± 12.3 (OP e) ± 4.3 {syst) ) MeV and f Bs = (232.5 ± 18.6 (OPE ) ± 2.4 (syst) ) MeV. 
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INTRODUCTION 



43 

The extraction of the ground-state decay constant within the method of QCD sum rules [l| is a complicated problem: 
First, one should derive a reliable operator product expansion (OPE) for the Borel-transformed (vulgo "Borelized") 
correlator II(t). The OPE in QCD is a double expansion: a perturbative expansion in powers of the strong coupling a s 
and an expansion in powers of the Borel parameter r in terms of condensates of increasing dimension. In practice, one 
has a truncated double series for which at most a few lowest-order terms can be calculated. This truncated series for the 
correlator depends explicitly on the renormalization scheme and scale (even in those cases where the full correlator does 
not), and the magnitude of the unknown higher-order terms crucially depends on the relevant choice of this scheme and 
\Q . scale. Therefore, controlling higher-order perturbative corrections poses a serious problem. 

CS| ' Since the pioneering work 0] , the correlator expressed in terms of the on-shell (or pole) heavy-quark mass has been 
. employed for the extraction of the decay constant. However, after the three-loop result for the correlator has appeared 
[H it became evident that the perturbative expansion in terms of the on-shell mass shows no signal of convergence: the 
LO, NLO, and NNLO terms all give comparable contributions to the decay constant. In contrast to this, a reshuffling of 
the perturbative expansion making use of the MS mass of the heavy quark leads to a clear hierarchy of the perturbative 
contributions to the decay constant [I]. 1 Following Q, we adopt the OPE of the relevant correlator in terms of the MS 
heavy-quark mass, denoted hereafter by uiq = rnQ(m,Q). Therefore, the on-shell mass does not appear in our analysis. 

Second, the knowledge of the correlator for only moderate values of r allows one to extract the characteristics of the 
bound state with some error which reflects the intrinsic uncertainty of the method of QCD sum-rules. Gaining control 
over this systematic uncertainty is a rather subtle problem, as it has been shown in 6] . Moreover, since higher multiloop 
perturbative calculations are becoming available and the knowledge of the fundamental QCD parameters is improving, 
the accuracy of the OPE of the relevant correlators is increasing. Therefore, the intrinsic systematic uncertainty of the 
QCD sum-rule method may become competitive with the decreased OPE uncertainties, as we shall show to be the case 
for the D-meson decay constant fjj. 

This work presents a detailed analysis of the decay constants of the D( s ) and mesons, with emphasis on acquiring 
control over all the uncertainties — of both OPE and intrinsic (i.e., systematic) origin — in these quantities. Recently, 
we formulated a novel algorithm for extracting bound-state parameters from OPEs for the correlators which opens the 
possibility to arrive at realistic error estimates for the extracted hadron parameters [7[ . The efficiency of our algorithm 
has been established in potential models: there the exact ground-state decay constants may be computed by solving the 



As rather unpleasant consequence, the decay constant extracted from the three-loop correlator in terms of the on-shell mass turns out to 
be considerably smaller than the estimate obtained from the three- loop correlator in terms of the running MS mass For the translation 
from one scheme to the other one employs the three-loop relation between the on-shell and the running quark masses. 
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Schrodinger equation. Moreover, it has been explicitly demonstrated that the extraction procedures of the ground-state 
parameters in QCD and in potential models are very close to each other quantitatively as soon as the (quark-hadron) 
duality is implemented in both theories in the same way 

This paper is organized as follows: Section [5] summarizes existing results on the OPEs for pseudoscalar heavy-light 
two-point functions and presents details of our algorithm for extracting the ground-state contribution to this correlator. 
Section [3] sketches our analysis of fn and Jd s , testing and proving the efficiency of our formalism. Section |4] provides 
the corresponding analysis of fs and fs a - Section [5] summarizes our conclusions. 



2. CORRELATOR AND SUM RULE 



We consider the correlator 

II(p 2 ) = i J d 4 xe^(0\T (j 5 (x)jt(0)) |0> (2.1) 
of two pseudoscalar heavy-light currents 

j 5 (x) = (rriQ + m)q(x)i-f 5 Q(x). (2.2) 
The Borel-transformed OPE series for this correlator is of the form 

oo 

II(r) = J dse~ ST ppcrt(s,/i) + n powcr (r, fi), (2.3) 

(mQ+m) 2 



where the perturbative spectral density may be found as a power series in the strong coupling a s : 

Ppcrt ( S , fi) = ,<»>(,) + «=^pW(.) + ( p^(s) + -... (2.4) 

TT \ TT J 

The Borel-transformed correlator (|2.3p does not depend on the renormalization scale fi; however, both the perturbative 
expansion truncated to some fixed order in a s and the truncated power corrections depend on pi. Moreover, the relative 
magnitudes of the lowest-order contributions strongly depend on the choice of the renormalization scheme and/or scale. 

A crucial issue for the reliability of a truncated perturbative expansion is the magnitude of the unknown higher-order 
corrections. One can play around with the choice of the renormalization scale in order to obtain some properties of the 
known terms of the perturbative expansion. For example, one may choose the scale p, by minimizing the highest-order 
known correction or by requiring some hierarchy of the known perturbative contributions. Unfortunately, even a clear 
hierarchy of several lowest-order perturbative corrections does not mean that the subsequent corrections are also small. 
Very often, the variation of the renormalization scale p in some range is used as an attempt to probe the magnitude of 
the unknown higher-order corrections. We shall pursue this strategy too, although there seems to be no rigorous way to 
estimate the size of these corrections without explicitly calculating them. 

In all expressions below, the quark masses toq and m and the strong coupling a s denote the respective MS running 
quantities at the scale p. Note that (s) depends on p implicitly through the quark masses, whereas all higher-order 
spectral densities p( n '(s), n > 1, depend on p implicitly through the quark masses and contain, in addition, explicitly 
/i-dependent logarithmic terms. Both perturbative spectral density and power corrections are given below for fi — mq. 

A. Perturbative spectral density 

The leading-order spectral density is well-known: 
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For the spectral density of order a s we make use of the first two terms of its expansion in small mass to derived in [J]: 2 
P W (s) = p^)+P^(s) + 0(m 2 ), 

p% {s) = -^C F {m Q + to) 2 s(1 - x) {(1 - x) [4L 2 (x) + 21n.xln(l - x) - (5 - 2x) ln(l - x)] 

+ (1 - 2.t)(3 - x) hxx + (17 - 33a;)/2} , 
p W(s) = —^c F (m Q +m) 2 m Q m{(l-x) [4L 2 (x) + 21na;ln(l - x) - 2(4 - x) ln(l - a;)] 

+2(3-5.x + x 2 )lnx + 2(7-9a;)}, (2.6) 
where x = to,q/s, N c — 3, and CV = (iV 2 — l)/2N c . The order-a 2 spectral density reads 

p (2) (s) = i? (2) ' s + A 1/9 (2) + A 2 p (2) +0(to). (2.7) 

Here, R^' s is the spectral function defined by Eq. (8) of Ref. |3|, which is provided by the authors through the publicly 
available program rvs.m. The authors of [1| have calculated the spectral function p pcr t(s) for the case m = in terms 
of the heavy-quark on-shell mass. Rewriting the O(l) and 0(a s ) spectral densities i?(°)' s and R^' s of @ in terms of 
the running mass generates the corrections Aip^ and A 2 p^ to the spectral density p^(s). The explicit expressions 
for these corrections are given by Eqs. (14) and (15) of Ref. and will not be reproduced here. 



B. Power corrections 



For the power corrections we also make use of the expression from Q: 

n power (r, p = m Q ) = (toq + to) e m Q T 

2C F a s ( \ m Q T \ 



(2.8) 



X { -m Q (qq) 



1 



1 



TO . 2 , TO 2 o n TOnT / 

2^ (1 + ^ T) + -^ T + ^ 1 



1 /a 



12 \ 7T 



-GG 



The parameter Toq describes the dimension-5 mixed quark-gluon condensate Q. It is worth noticing that the radiative 
corrections to the condensates increase rather fast with the Borcl parameter r. 

In summary, we make use of the expressions for II(t) from Q with minor modifications: 

(i) We adopt the "natural" threshold (toq+to) 2 in the spectral representation for II(t) and therefore encounter only 
running MS masses in our formulas. 

(ii) For pW(s) we use the exact expression without performing an expansion in powers of to, and for p^\s) we do 
not include terms of order to 2 and higher. 

The OPE parameters required for our analysis are 

m(2 GeV) = (3.5 ± 0.5) MeV, m s (2 GeV) = (100 ± 10) MeV, 

, 3 ( S " S )(2GeV) 



(qq)(2 GeV) = -((267 ±17) MeV) 
m 2 = (0.8 ± 0.2) GeV 2 , a s (M z ) = 0.1176 ± 0.0020. 



= 0.8 ±0.3, (— GG) = (0.024 ±0.012) GeV 4 
{qq}(2 GeV) \ 7r 



(2.9) 



We perform the calculations for two sets of c- and 6-quark masses toq = toq (toq): the values from PDG Q 

m c = (l.27±g;SX) GeV, m b = (4.2±°;^) GeV, (2.10) 
and the very accurate values reported recently in fioj 

(1.279 ±0.013) GeV, mb — (4.163 ± 0.016) GeV. (2.11) 



rn c 



2 Note that the 0(a a m 2 ) corrections to II(t) cannot be obtained by expanding the spectral density p' 1 ) in powers of m: Starting from the 
order m 3 , the functions p^ 3 etc. contain poles of increasing orders at s = rrig (see Eqs. (C3) and (C4) of Ref. Q). Therefore, after the 
s-integration all these terms yield contributions of the same order 0(a s m 2 ) to the spectral function. For the same reason, the expansion 
of (s) in powers of m does not allow one to get the terms of order m 3 in n(r). We therefore use the exact expression for the spectral 
density p(°) instead of expanding it in powers of m. 
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C. Sum rule 

The correlator (|2.ip may be evaluated by inserting a complete set of hadronic intermediate states: 

n(r) = ]qM^ m Q t + contributions of higher states, (2.12) 

where Jq is the decay constant of the Pq meson, defined according to 

(m Q + m)(0\ui l5 Q\P Q ) = f Q M 2 Q . (2.13) 

For large values of r, the contributions of the excited states decrease faster than the ground-state contribution and the 
correlator II(t) is dominated by the ground state. Unfortunately, the truncated OPE does not allow us to calculate the 
correlator at sufficiently large r, such that the excited states give a sizable contribution to n(r). 

According to the quark-hadron duality assumption, the contributions of excited states and continuum are described 
by the QCD perturbative contribution above an effective continuum threshold s c g- This leads to the following relation: 

Scff(-r) 

/qMje- M « r -n dual (r,s cff (r)) = J ds e~ ST p port (s, fi) + n powcr (r, //). (2.14) 

(mQ+m) 2 

In the region near the physical continuum threshold at s — (M\/ Q +m v ) 2 , Vq being the lightest vector meson containing 
a quark Q, the QCD perturbative spectral density and the hadron spectral density are rather different. Consequently, 
the effective continuum threshold as defined by (|2.14j) turns out to be necessarily a function of the Borel parameter r. 
We introduce the dual invariant mass -Mduai and the dual decay constant /dual by the definitions 

-^-logn dual (r, Scff (r)), (2.15) 

CZT 

M- 4 e M « T n dual (r, Scff (r)). (2.16) 

Notice that the deviation of the dual mass from the actual ground-state mass provides an indication of the excited-state 
contributions picked up by the dual correlator. 

In order to determine the decay constant /q of the Pq meson from the OPE we must execute the following two steps. 

1. Borel window 

First, we have to fix the working r-window where, on the one hand, the OPE yields a sufficiently accurate description 
of the exact correlator (that is, the higher-order radiative and power corrections are small) and, on the other hand, the 
ground state gives a "sizable" contribution to the correlator. Since the radiative corrections to the condensates increase 
rather fast with r, it is preferable to stay at the lowest possible values of r. We shall therefore fix the Borel window by 
the following criteria: (a) In the window, the power corrections IipowcrM should not exceed 30% of the cut perturbative 
correlator n pcr t(r, so); this gives the upper boundary of the r-window. The ground-state contribution to the correlator 
at such r values comprises about 50% of the correlator, (b) The lower boundary of the r-window is defined by requiring 
the ground-state contribution not to drop below 10%. In quantum physics, our algorithm was shown to provide a good, 
even excellent extraction of the ground-state decay constant for Borel windows determined by these requirements [6j,[7|]. 

2. Effective continuum threshold 

Second, we must formulate our criterion for the determination of s e g(T). We consider an algorithm for the extraction 
of Jq which takes advantage of the knowledge of the PQ-meson mass Mq. Our algorithm, constructed in previous works 
and, in quantum-theoretical potential models, proven to work well for various correlators, is rather simple: We consider 
a set of r-dependent Ansatze for the effective continuum threshold, for simplicity assumed to be all of polynomial form: 

n 

4ff ) (r) = E4" ) ^- ( 2 - 17 ) 

We determine the parameters on the r.h.s. of (|2.17p as follows: We calculate the dual mass squared according to (|2.15[) 
for the r-dependent s c g of Eq. (|2.17l) . We then compute Af| ual (r) at several values of r = r, (i = 1, . . . , N, where N can 



M| ual (r) = 

/dual( r ) — 
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be taken arbitrary large) chosen uniformly over the Borel window. Finally, we minimize the squared difference between 
M| ual and the known value Mq: 

i N 

This pins down the parameters of the effective continuum threshold. As soon as the latter is fixed, it is straightforward 
to calculate the decay constant /q. 

According to our recent findings, allowing for a r-dependence of the effective threshold leads to visible improvements 
compared with the traditional assumption of a r-independent quantity: The former yields a much better stability of the 
dual mass calculated from the dual correlator and allows one to work at smaller values of r, where the impact of power 
corrections is reduced. 



3. Uncertainties of the extracted decay constant 

The above discussion implies that the extracted result for the decay constant is sensitive both to the precise values of 
the OPE parameters and to the particular prescription for fixing the effective continuum threshold. The corresponding 
uncertainties of the predicted decay constant are labeled as its OPE-related error and its systematic error, respectively: 

f . OPE-related error. The OPE-related uncertainty is estimated as follows: We perform a bootstrap analysis [ll[ by 
allowing the OPE parameters to vary over the ranges quoted in Eqs. (|2.9p . using f 000 bootstrap events. Gaussian 
distributions for all parameters but [i are employed. For fi we assume a uniform distribution in the corresponding 
range, which we choose to be f < [i (GeV) < 3 for charmed mesons and 2 < fi (GeV) < 8 for beauty mesons. The 
resulting distribution of the decay constant turns out to be close to a Gaussian shape. The quoted OPE-related 
error is therefore the Gaussian error. 

2. Systematic error. The systematic uncertainty of some hadron parameter obtained by the sum-rule method (i.e., 
the error related to the intrinsically limited accuracy of this approach) represents the perhaps most subtle point 
in all applications of this method. So far no way to provide a rigorous — in the mathematical sense — systematic 
error has been devised. However, a realistic estimate of the corresponding error may be found: As prompted by 
detailed comparisons of the extraction of the decay constant in QCD and in potential models [§[ , the band of fp 
values spanned by the linear, quadratic, and cubic Ansatze for the effective continuum threshold contains the true 
value of the decay constant. Trusting in these findings, the half-width of this band is interpreted as the systematic 
error of the decay constant. Presently, we do not envisage any other possibility to provide more reliable estimates 
for the systematic error. 



3. DECAY CONSTANTS OF THE D AND D 3 MESONS 



A. Decay constant of the D meson 

The Borel window for the charmed meson is chosen according to the criteria discussed above: r = (O.f — 0.5) GeV -2 . 
Figure [T] illustrates the application of our prescription of obtaining the effective continuum threshold and extracting the 
corresponding fp>. We would like to point out that r-dependent effective thresholds lead to a much better reproduction 
of the meson mass in the window than a constant one (Fig.QJ,). This signals that the dual correlators corresponding to 
the r-dependent thresholds are less contaminated by the excited states. 

The dependence of the extracted value of the D-meson decay constant fu on the c-quark mass m c = m c (rn c ) and the 
quark condensate (qq) = (qq) (2 GeV) may be parameterized in the form 



20<,2 - ! 1279 M + 4 r Kg " g)|1 ^7°- 267 GeV 1 ± *. W, . 



( \mv 

O.f GeV J ' ~\ O.Of GeV 



MeV. (3.1) 



This relation describes the band of values delimited by the two dotted lines in Fig.[Tb, which include the results derived 
with the linear, quadratic, and cubic Ansatze for the effective continuum threshold. Figure^ displays the results of the 
bootstrap analysis of the OPE uncertainties. The distribution has a Gaussian shape, and therefore the corresponding 
OPE uncertainty is the Gaussian error. Adding the width of the band provided by the r-dependent n = 1, 2, 3 Ansatze 
for the effective continuum threshold as the (intrinsic) systematic error of our approach, we obtain the following result: 

f D = (206.2 ± 7.3 ( ope) ± 5.1 (syst) ) MeV. (3.2) 
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Fig. 1: Dual mass (a) and dual decay constant (b) of the D meson extracted by adopting different Ansatze (|2.17[) for the effective 
continuum threshold s e ff (r) and fixing these thresholds according to (|2.f 8[) . The results for m c = m c (m c ) = 1.279 GeV, /i = m c , 
and central values of the other relevant parameters are presented, (c) Dual decay constant of the D meson vs. m c for /i = m c and 
central values of the other OPE parameters. The index n = 0, 1, 2, 3 denotes the power of our polynomial Ansatz for the effective 
continuum threshold in (12.171) . 
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Fig. 2: Distribution of fo as obtained by a bootstrap analysis of the OPE uncertainties. For all OPE parameters but fi Gaussian 
distributions with corresponding errors quoted in (|2.9p are employed. A variation of m c in the interval (I2.10p (a) and (|2.1ip (b) 
is allowed. For we assume a uniform distribution in the range 1 GeV < fj. < 3 GeV. (c) Summary of our results for fn ■ For 
comparison, the lattice results are shown for two dynamical light flavors (Nf = 2) [l^, El an d three dynamical flavors (TV/ = 3) 
[3 EH- The triangle represents the experimental result from PDG For the r-dependent QCD-SR result the depicted error is 
the sum of the OPE and systematic uncertainties given in (|3.2p . added in quadrature. 



Wc have considered for m c the two ranges in (|2 . 1 0[) and (|2 . 1 1 1) . The OPE- related error is practically the same for both 
ranges (see Figs. [2Kb), so the main source of the OPE uncertainty in the extracted fo comes from the OPE parameters 
other than m c (mainly, the renormalization scale and the quark condensate) . 

Notice that the bootstrap procedure for a r-independent effective threshold gives a substantially lower fu range, viz., 
/d(« = 0) = (181.3±7.4(ope)) MeV, which deviates from our r-dependent result (|3.2p by almost three times the OPE 
uncertainty. Moreover, as we have already shown in our previous works 0], making use of merely a constant Ansatz for 
the effective continuum threshold docs not allow one to probe at all the intrinsic systematic error of the QCD sum rule. 
From (|3.2p the latter turns out to be of the same order as the OPE uncertainty. 

The r-dependent threshold leads to a clearly discernible effect and brings the results from QCD sum rules into perfect 
agreement with recent lattice results and the data (Fig. (5}d) . A perfect agreement of our result with the lattice ones and 
with experiment provides a further confirmation of the reliability of our procedure. 



B. Decay constant of the D 3 meson 



The corresponding r- window is r — (0.1 — 0.6) GeV 2 . Figure[3]provides the details of our extraction procedure. Our 
results for fu a may be represented in the form 



7 




Fig. 3: Same as Fig. \T\ but for the Ds meson. 
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This relation describes the band of fo 3 values as a function of m c = m c (rn c ) indicated by the two dotted lines in Fig. [3};, 
as well as the dependence on the quark condensate (ss) = (ss)(2 GeV). Performing the bootstrap analysis of the OPE 
uncertainties, we obtain the following estimate: 

f Ds = (245.3 ± 15.7 (0 pe) ± 4.5 (syst) ) MeV. (3.4) 

As in the case of /d , a constant threshold yields a substantially lower fo value: fD s {n = 0) = (218.8±16.1(qpe)) MeV. 



C fnjfo 

For the ratio of the D( s )-meson decay constants, we find 

IdJId = 1-193 ± 0.025 (O PE) ± 0.007 (syst) , 



(3.5) 



to be compared with the PDG average, fp / fp = 1.25±0.06 and the recent lattice results fpj fp = 1.24±0.03 [12| 
at N f = 2, and f D J Id = 1.164±0.011 [3] and f D Jf D = 1.20±0.02 [U] at N f = 3. The erro/in comes mainly 
from the uncertainties in the quark condensates: (ss)/(qq) — 0.8 ± 0.3. 
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Fig. 4: (a) Distribution of /d s as obtained by the bootstrap analysis of the OPE uncertainties. For all our OPE parameters but 
fi Gaussian distributions with corresponding errors given in (12. 9[) are adopted. A variation of m c in the interval (|2.11|l is allowed. 
For n we assume a uniform distribution in the range 1 GeV < n < 3 GeV. (b) Summary of our results for Jds- Lattice results are 
depicted for N f = 2 [H, [T3] and N f = 3 GJ] . The experimental results are from PDG and the CLEO collaboration Hi]. 
For the r-dependent QCD-SR result the error shown is the sum of the OPE and systematic uncertainties given in l|3.4fl . added in 
quadrature. 



8 



DECAY CONSTANTS OF THE B AND B s MESONS 



We set the Borel window as r = (0.05 — 0.18) GeV 2 . Note that the radiative corrections to the condensates increase 
rather fast with r, so it is preferable to stay at lower values of t. 



A. Decay constant of the B meson 



Figure [5] shows the application of our prescription to the extraction of fg. The correlator II(t), which has dimension 
six, is extremely sensitive to the precise value of mj. The dependence of the extracted value of the decay constant fg on 
the 6-quark mass m b = m b (m b ) and the condensate (qq) = (qq)(2GeV) may be parameterized in the form 



m b , {qq}) 



( I (qq) 1 1/3 - 0.267 GeV 



0.01 GeV 



±4 



(syst) 



MeV. (4.1) 



The above relation describes the band of values delimited by the two dotted lines in Fig.JSJ: which, as before, include the 
results obtained with the linear, quadratic, and cubic Ansatze for the effective continuum threshold. In addition, it also 
encodes the dependence on the value of the quark condensate. 




Fig. 5: Dual mass (a) and dual decay constant (b) of the B meson obtained by using different Ansatze for the effective continuum 
threshold s e ff(r) (|2.17(l and fixing the coefficients according to (|2.18|l . The results for rrib = mb(rfib) = 4.245 GeV, /i = nib, and 
central values of all other relevant parameters are presented, (c) Dual decay constant of the B meson vs. mi, for [i = mj, and 
central values of all other OPE parameters. The index n = 0, 1, 2, 3 denotes the power of the polynomial Ansatz for the effective 
continuum threshold in (|2.17[) . 



We now perform a bootstrap analysis of fg combining all OPE uncertainties. We assume Gaussian distributions for 
the OPE parameters (quark masses, condensates) with corresponding errors. The renormalization scale /i is assumed to 
be uniformly distributed in the interval 2 < [i (GeV) < 8. 

Because of the high sensitivity of the correlator to the fo-quark mass, the sum-rule estimate for fg strongly depends on 
the range of mb used. For the PDG range m b — (4.2+q'qy) GeV Q, the OPE uncertainties are very large; therefore, no 
reasonable estimate of fg may be obtained (see Fig. [6]). 

Adopting the recently reported precise range m b = (4.163±0.016) GeV from Qjl leads to a rather accurate estimate: 
fg = (225.6±11.3(opE)±2.2( syst )) MeV. However, we observe (see Fig. [7]) some tension between this value and recent 
lattice calculations of fg, which yield, on average, /^(lattice) = (193 ± 130) MeV. Requiring the sum-rule estimate to 
match the lattice average leads to a rather accurate determination of m b : 

m b (m b ) = (4.245 ± 0.025) GeV. (4.2) 

This, however, differs considerably from the range found in fiol ] as well as from the recent finding rn b (m b ) = 4.164 ± 
0.023 GeV [2Q], obtained from a (perturbative) QCD analysis similar to the one used in [l(| but applied to the 
moments of heavy-quark current-current correlators calculated in lattice QCD with Nf = 3 . The value (14.21) is 
in good agreement with the lattice determinations m b (m b ) = 4.26 ± 0.03 s ^ a ^ ± 0.09 S y S ^ GeV [2l| and m b (m b ) = 
4.25 ± 0.02 s t a t ± 0.11 S y S t GeV [Hf, as well as with the recent preliminary result of the Alpha Collaboration [2j|, all 
of them obtained using HQET on the lattice with Nf = 2. 

The B-meson decay constant and its uncertainties corresponding to the 6-quark mass (|4.2|) are 

fg = (193.4 ± 12.3 (0 pe) ± 4.3 (syst) ) MeV. (4.3) 

Finally, let us mention that, for the range (14.21) . a constant effective threshold gives f B (n — 0) = (184 ± 13(qpe)) MeV. 
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Fig. 6: Distribution of /b for three ranges of mi,: the range (|2.10[1 (a), the range (12. lip (b), and the range specified by (|4.2[) (c). 
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B. Decay constant of the B a meson 

Figure [S] depicts the application of our procedure to the extraction of fs s ■ The dependence of the extracted value of 
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Fig. 8: Extraction of /s s . Line identification as in Fig. [5] 

the decay constant fg a on the 6-quark mass and the strange-quark condensate is given by the relation 



m b [GeV] 



0.1 GeV J \ 0.01 GeV 



(syst) 



MeV. (4.4) 
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This formula describes the band of values indicated by two dotted lines in Fig. [8): and, in addition, gives the dependence 
on the value of the quark condensate at renormalization scale fj, = 2 GeV. 

We now perform a bootstrap analysis of Jb s combining all OPE uncertainties. We assume Gaussian distributions of 
the OPE parameters (quark masses, condensates) with corresponding errors. The renormalization scale /i is assumed to 
be uniformly distributed in the interval 2 < /i (GeV) < 8. Making use of the range mj = (4.163 ± 0.016) GeV [10| one 
has f Bs = (262.0 ± 18.1 ( ope) ± 2.9 (syst) ) MeV, while the range g2J), m b = (4.245 ± 0.025) GeV, leads to 



J B . = (232.5 ± 18.6 (OP E) ± 2.4 (syst) ) MeV. 
Figure [SJa compares our results with recent lattice determinations. 
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Fig. 9: (a) Distribution of /s s for mj = (4.245 ± 0.025) GeV [sec (|4.2|| ]. (b) Our results, corresponding to the m(,-ranges (|2.11[) 
and UJJ, vs. lattice-QCD results for N f = 2 [H, EE3] and N f = 3 |lElIl- 



C fsjfs 

For the ratio of the decay constants of beauty mesons, we find 

SbJSb = 1-203 ± 0.020 ( ope) ± 0.007 (syst) , 



(4.6) 



which has to be compared with the lattice results ! b J!b = 1.27±0.05 at N f = 2 [Tg and JbJIb = 1.226±0.026 (3 
and JbJ 'Jb = 1.245±0.043 [19( at Nf = 3. Similar to the charmed ratio (|3.5p . the error in (|4.6[) comes mainly from the 
uncertainty in the ratio of quark condensates: (ss) / (qq) — 0.8 ± 0.3. 



The double ratio of the decay constants — the Grinstein ratio R\ [2 
We obtain the particularly accurate value 



D. Double ratio (f Ba /fB)/(f D Jfo) 

- is practically free from OPE uncertainties. 



IbJ ft 



1 = 0.013 ±0.011 



(syst); 



Jd 3 I Id 

which is consistent with the lattice determination (f B J / B )/(/rj 8 / f D ) = 0.018 ± 0.006 ± 0.010 [H]. 



(4.7) 
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5. SUMMARY AND CONCLUSIONS 



In summary, we performed a detailed analysis of the extraction of the decay constants of pseudoscalar heavy mesons 
from the correlator of pseudoscalar currents. Particular emphasis was laid on the investigation of the uncertainties in 
the extracted values of the decay constants: namely, on the OPE uncertainty related to the not precisely known QCD 
parameters, and on the intrinsic uncertainty of the method related to the limited accuracy of the extraction procedure. 
According to our recent findings, the accuracy of the sum-rule estimates may be considerably improved and the intrinsic 
uncertainties in hadron parameters may be probed by studying systematically the Borel-parameter dependence of the 
effective continuum thresholds; the parameters of these effective thresholds may be fixed by minimizing the deviation of 
the dual mass from the known meson mass in the Borel window. In the present work, this strategy has been applied to 
the decay constants of heavy mesons. Our main results are as follows: 

(i) We obtain the following estimates for the decay constants of the charmed D and D s mesons: 

f D = (206.2 ±7.3 (0 PE)±5.1 (syst) )MeV, (5.1) 
f Da = (245.3 ±15.7 (0 P E )±4.5 (syst) )MeV. (5.2) 

We would like to point out that we provide both the OPE uncertainties and the intrinsic (systematic) uncertainty of the 
method of sum rules related to the limited accuracy of the extraction procedure. In the case of /d, the latter turns out 
to be of the same order of magnitude as the OPE uncertainty. Noteworthy, assuming r-independence of the effective 
continuum threshold leads to the substantially lower decay-constant range fo{n = 0) = (181.3±7.4(qpe)) MeV, which 
differs by almost three times the OPE uncertainty from our result (|5.ip found from a r-dependent effective threshold. 
The ratio of the charmed-meson decay constants (|5.2p and (|5.1[) is 

f D Jf D = 1-193 ± 0.025 ( ope) ± 0.007 (syst) . (5.3) 

(ii) The decay constants of the B and B s mesons are very sensitive to the precise value of nib- Using the PDG range 

of nib does not allow us to obtain a reasonable estimate. For the very narrow range mb(mb) — (4.163±0.0016) GeV [lOj . 

our analysis gives f B = (225.6 ± 11.3 (0 p E ) ± 2.2 (syst) ) MeV and f Bs = (262.0 ± 18. l (0 p E ) ± 2.9 (syst ))MeV. We observe 

some tension between the above sum- rule result for fg and the average of recent lattice calculations (13l Ii~7l - fl9| . namely, 
^lattice) = lg3 ± 13 MeV ^ 

We emphasize that the observed strong sensitivity of f B to the precise value of nib provides an interesting alternative 
way of obtaining nib from the analysis of the decay constant: Using the lattice average for f B as input yields the rather 
accurate estimate for the 6-quark mass 

m b (m b ) = (4.245 ± 0.025) GeV. (5.4) 

This new range of m b corresponds to 

f B = (193.4 ± 12.3 ( ope) ± 4.3 (syst) ) MeV (5.5) 

and yields 

f Bs = (232.5 ± 18.6 (0 pe) ± 2.4 (syst) ) MeV. (5.6) 
For the ratio of the decay constants (15.61) and (I5.5[) we get 

SbJSb = 1-203 ± 0.020 ( ope) ± 0.007 (syst) , (5.7) 

(iii) The double (Grinstein) ratio of the decay constants, 

m^--l = 0.013 ±0.011 (syst)) (5.8) 
JDs I Id 

is practically free from OPE uncertainties and, consequently, may be predicted with rather high accuracy. 
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